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We perform a field-theoretical computation of hadron pro- 
duction in large systems at the QCD confinement phase tran- 
sition associated with restoration of the Z{3) global sym- 
metry. This occurs from the decay of a condensate for the 
Polyakov loop. From the effective potential for the Polyakov 
loop, its mass just below the confinement temperature Tc is 
in between the vacuum masses of the pion and that of the 
kaon. Therefore, due to phase-space restrictions the number 
of produced kaons is roughly an order of magnitude smaller 
than that of produced pious, in agreement with recent results 
from collisions of gold ions at the BNL-RHIC. From its mass, 
we estimate that the Polyakov loop condensate is character- 
ized by a (spatial) correlation scale of 1/me — 1/2 fm. For 
systems of deconfined matter of about that size, the free en- 
ergy may not be dominated by a condensate for the Polyakov 
loop, and so the process of hadronization may be qualitatively 
different as compared to large systems. In that vein, exper- 
imental data on hadron abundance ratios, for example K/n, 
in high-multiplicity pp events at high energies should be very 
interesting. 



High-energy inelastic processes provide the unique op- 
portunity to recreate perhaps the high energy density 
state of QCD matter that presumably prevailed during 
the first microseconds of the evolution of the early uni- 
verse: the so-called Quark-Gluon Plasma (QGP). Sim- 
ply speaking, the QGP is the state of matter for which 
quarks and gluons are deconfined and chiral symmetry is 
restored. Even if produced, the QGP is a transient state 
of matter which can not be observed directly. Experi- 
mentally, one observes hadrons which are produced from 
the decay of the deconfined QGP state. The hadron pro- 
duction process itself is commonly called hadronization. 

One might be able to gain some insight into the 
hadronization process by studying the relative popula- 
tions of various hadronic states emerging from a high- 
energy inelastic process 0. In fact, one of the earliest 
proposed signatures for a transient QGP state is the en- 
hanced production of strangeness, particularly the rela- 
tive abundance of kaons to pions (the "i^/TT ratio") §. 
More recently, it has been discussed that a large variety of 
ratios of hadronic multiplicities in high-energy heavy-ion 
collisions resemble those of a hadron gas at a tempera- 
ture Th close to Tc, the confinement/deconfinement tran- 
sition temperature (^-|^ . This has been interpreted such 
that the QGP state hadronizes statistically, in that all 



hadronic states are populated according to their weight in 
the partition function of a hadron gas [§,0J|] ■ This inter- 
pretation was partly based also on the success of the sta- 
tistical hadronization model to describe relative hadron 
abundances in high-energy inelastic e'^e~ and pp, or pp 
reactions In such reactions, one does not expect that 
a thermally and chemically equilibrated QGP is formed 
prior to hadronization. Nevertheless, the hadronic final 
states appear to be populated approximately according 
to a statistical distribution. From e+e~ to central colli- 
sions of large nuclei, the hadronization temperature ob- 
tained from fitting the relative multiplicities was found 
to be more or less the same, and that it is also indepen- 
dent of collision energy at high energies This was 
interpreted as a "limiting temperature" for hadron 
formation in inelastic reactions, which is consistent with 
the prediction of a deconfinement temperature Tc from 
lattice QCD above which hadronic states presumably can 
not freeze out as asymptotic states. 

Here, we perform a dynamical computation of hadron 
production from the decay of a deconfined state with a 
spontaneously broken global Z{3) symmetry without as- 
suming the formation of an equilibrated hadron gas a 
priori. For simplicity, we restrict ourselves to the pseu- 
doscalar meson octet, i.e. the pions, the kaons, and the 
eta. At the end of the paper, we also speculate about 
possible differences in the "statistical hadronization" in 
small versus large systems. 

We begin by briefly reviewing the model for the QGP 
suggested recently by Pisarski which we shall adopt 
here to compute hadron production. The basic postulate 
of the so-called "Polyakov loop model" is that the free 
energy of the deconfined state of QCD is dominated by 
a condensate of gauge invariant Polyakov loops 
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It is based on the observation from SU{?>) lattice 
gauge theory that the Polyakov loop becomes light near 
Tc 0,0. In other words, the effective potential for £ 
must be rather flat for T ~ Tc- On the lattice, the mass 
of £ can be measured from the two-point function with 
space-like separation. These two point functions in the 
Polyakov loop model are very different from those of or- 
dinary perturbation theory . The condensate is then 
characterized by a length scale ^ = l/me{T). 

In order to study hadronization, we couple the 
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Polyakov loop to the SU{3) x S'L''(3) linear sigma model 
and make qualitative predictions for the resulting relative 
multiplicities of kaons and etas to pions. The results of 
this study have implications for understanding the flavor 
composition in the final state of high-energy inelastic re- 
actions. For simplicity, we neglect the effects of nonzero 
baryon charge, isospin charge, electric charge, etc. in the 
present analysis. This should be a reasonable first ap- 
proximation for particles produced in the central region 
of very high energy (and high multiplicity) inelastic pro- 
cesses. In other cases, one may extrapolate the measured 
hadronic multiplicities to vanishing charges |^ , for exam- 
ple, so as to allow for a better comparison of e+e^, pp, 
and nucleus- nucleus collisions. 

Our model has static, thermodynamic properties which 
are consistent with lattice data for the pressure, the en- 
ergy density, and the two-point function (i.e. the mass) of 
the Polyakov loop. With the simplest possible Z(3) sym- 
metric potential described below, one can not only fit the 
above observables but actually relate their behavior, for 
example the pressure above Tc to the electric screening 
mass [|2|. One also obtains testable predictions for the 
ratio of the masses of the real and imaginary parts of the 
Polyakov loop . 
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FIG. 1. The energy density and the pressure for SU(3) 
pure-gauge theory, and our fit (see text). 

For three colors, we take the potential for £ to be 
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which is invariant under global Z{3) transformations 
1^,^. The form of V{i) is dictated by symmetry princi- 
ples [Hn). The center of SU{3), Z(3), is broken at high 
T where the global minimum of V{i), f-o(T) — > I, and 
restored at low T where ^o(T) — > 0. 

We assume that quarks are not important for under- 
standing the form of the effective potential for QCD. This 
is motivated by the results of lattice simulations for which 
P/Pideai as a function of T/Tc is remarkably insensitive 



to the number of flavors, Nf = 0, 2, 3 [Q. Thus, we ne- 
glect terms linear in £, which exist in QCD with quarks. 
Such terms would change the weakly first-order transi- 
tion (for Nf = 0) which is predicted by eq. (||) into a 
crossover. Nevertheless, the fact that the free energy is 
small just below Tc means that their numerical effects 
must be small. We emphasize that the results to be de- 
scribed below are not driven by the order of the phase 
transition and would not change in any significant way 
by the inclusion of such small linear terms. All that mat- 
ters is that the £ field becomes rather light at Tc, i.e. that 
the transition is nearly second-order. 

The coefficients b2{T), b^, and 64 can be chosen to 
reproduce lattice data for the pressure and energy den- 
sity of the pure glue theory for T > jl^. For ex- 
ample, one may chose 64 ^ 0.6061, 63 = 2.0, and 
b2{T) = (l-l.II/x)(l-H0.265/a;)2(l-h0.300/a:)3- 0.487 
||l3| , where x = T jT^- For this parameterization, one ob- 
tains the pressure and energy density e = Tdp/dT — p as 
shown in Fig. |^. The above parameterization does not yet 
incorporate the normalization condition that ^o(T) 1 
for T 00. Rather, the expectation value for £ at high 
temperature approaches r = 63/2-1-1/63 -I- 462(oo)/2. We 
thus rescale the field and the coefficients as ^ — > £/r, 
62(T) ^ 62(T)/r2, 63 ^ 63/r, and 64 ^ b^r'^. This 
ensures the proper normalization for the Polyakov loop. 
The rescaling does not, of course, change the pressure 
and the energy density. 




FIG. 2. The potential for the real part of £ at three different 
temperatures. 

For Nf — 3, we rescale 64 by a factor of (47 + ^)/f6, 
corresponding to the increase in the number of degrees 
of freedom. In the spirit of mean field theory, the same 
constant values of 63 and 64 are employed for T < Tc- 
In this domain, 62 (T) can be set using the string ten- 
sion cr(r) from the lattice g^/2Nc d^V/d£^ = 
-'g^b2{T)biT^/2Nc = <J^{T). The factor in front of 
V" arises from the normalization of the kinetic term 
for see below. Taking Nc/g^ = 1, this leads to 
62(T) « -0.66ct2(t)/T4. 
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The potential in eq. (g) , which follows from the above- 
mentioned fit to the lattice data, changes extremely 
rapidly below Tc see Fig. ||. This rapid change in 
the effective potential with temperature has important 
implications for the dynamical evolution of the confine- 
ment transition p^], and for the relative hadron mul- 
tiplicities from the decay of the condensate for ^. The 
hadron production temperature Th can not vary much 
between events because it is impossible for the expecta- 
tion value for t to stay large below Tc; rather, the field is 
in an unstable configuration and must roll down towards 
the minimum of the potential. (In turn, if the effective 
potential was slowly varying about Tc, one would expect 
a broad range of hadronization temperatures.) Also, Th 
must be very close to Tc, which in turn implies that the 
K/tt ratio is well below unity. As we shall see below, if 
we chose Th by hand to be some 10% or more below Tc 
(disregarding the fact that from Fig. || this scenario ap- 
pears unphysical), then phase-space for kaon production 
opens up and leads to essentially equal number of kaons 
and pions in the final state, contrary to experimental re- 
sults. Thus, the question as to whether it is some feature 
of QCD that the hadron production temperature is so 
close to the confinement temperature Tc is answered 
by the behavior of the potential for the Polyakov loop 
which we deduce from the lattice data. 

To complete the effective theory, we add a kinetic term 
for t and couple it to the octet of pseudo-Goldstone exci- 
tations. The coefficient of the kinetic term for £ is set by 
the fluctuations of the SU (3) Wilson lines in space. Per- 
turbative corrections to that coefficient are small [ p^ . 

The 



We assume here a Lorentz invariant form |13 
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coupHng of i to the SUr{3) x 5[/l(3) hnear sigma model 
is given by 

C = C^ + ^\dJ\^T^-V(i)-h^\e\'^T^TT i^'i''') (3) 

g 

where is the Lagrangian for the linear sigma model 
given below and $ is the chiral field. The coupling ~ 
20 between the chiral field and £ is chosen to reproduce 
77i7r (T) . Through this interaction term, the Polyakov 
loop condensate drives the chiral phase transition. 

The Lagrangian of the SU{3)ii x SU{3)l hnear sigma 
model is given by ||2^ 

= Tr (9^$taM$ _ ^2 ^t$) _ Xi [Tr (4) 

- AaTr + c [Det ($) + Det 

-I- Tr [H{<P + $t)] . 

$ is a complex 3x3 matrix parameterizing the scalar 
and pseudoscalar meson nonets. 
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where Aq are the Gell-Mann matrices with Aq = y 1 
The (Ja fields are members of the scalar nonet and the 
TTa fields are members of the pseudoscalar nonet. Here, 
tt"^ = (tti ± 11x2)1^/2 and tt" = 7r3 are the charged and 
neutral pions, respectively. = (7r4 ± iiT^)/^/2, = 
{ttq + inj)/^/2, and K'^ = {ttq — are the kaons. 

The 77 and the 77' mesons are admixtures of ttq and TTg. 
The classification of the scalar nonet is not important for 
our purposes as discussed below. 

The 3x3 matrix H is chosen to reproduce the mass 
splitting of the strange and nonstrange quarks with H — 
Ha\a/'2' and where Ha are nine external fields. The deter- 
minant terms reproduce the effects of the U{\)a anomaly 
in the QCD vacuum |^ by preserving the S'[/(3)fl x 
SU{i)L = SU{i)v X SU{i)A symmetry while explicitly 
breaking the U{1)a symmetry. Symmetry breaking gives 
the <i> field a vacuum expectation value, (<i>) = Ua^aj'^- 
The fJo and a% are admixtures of the strange and non- 
strange quark-antiquark condensates. From the PCAC 
relations, they are given by ctq = (/,r + j — 
130 MeV and ag = 2 (/^ - /k) /V3 = -23.8 MeV. Here, 
for simplicity we freeze tro, cts at their vacuum values, 
i.e. we assume that dynamical relaxation of the classical 
condensates does not contribute much to particle pro- 
duction or to the tadpole self-energies of the tt, r\ 
quantum fields. From the numerical values of ctq, 0% it is 
clear that the energy density of the chiral condensate is 
too low to contribute much to the total hadron multiplic- 
ities. The free energy of the deconfined state is largely 
due to the potential (||) for the Polyakov loop I which 
dominates hadron production at Tc IT^ . 

The parameterization of the coupling constants is given 
in detail in Ref. and is chosen to reproduce the masses 
and mixing angles of the mesons at zero temperature. 
Here, we take = (342 MeV)^, Ai = 1.4 and A2 = 46. 
The explicit symmetry breaking terms are given by = 
(286 MeV) ^ iJg = -(311 MeV) ^ and c = 4808 MeV. 

To simplify the solution of the equations of motion for 
t and the mesonic excitations, we do not consider the 
excitations of the scalar nonet or the 77', all of which 
have masses either near or above 1 GeV. We perform the 
standard decomposition of the meson quantum fields in 
terms of creation and annihilation operators times adia- 
batic mode functions 7rJ(t). Here, fc is the mode number. 
To make the interaction quadratic in the meson fields, 
we employ a mean field factorization which respects the 
fiavor-S'C/ (3) symmetry when fluctuations are large, 
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With these approximations, the interaction part becomes 



-/i^T^I^I^I 7r^ + X^ + 77^ 



(7) 
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where we introduced the short-hand notation 



(8) 



for the sum of the fluctuations. The prefactor of 
each term reflects the internal isospin degrees of free- 
dom. Assuming that the up and down quarks are mass- 
degenerate, the equations of motion for the mesonic mode 
functions in Minkowski space read 
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(10) 



m2-Ai((02)+a2+a|) 
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The equations of motion are now N^. (coupled) equa- 
tions for the mode functions, where is the number 
of modes. In the actual computations we propagated 
Nk = 1000 modes in parallel, taking a mode-spacing of 
Afc = /tt/IOO. We checked that the observables, like the 
total number of produced pions and kaons, were stable 
under reasonable variations of Nk- When discussing rel- 
ative abundances of tt, K, 77 below, we always refer to 
the sum of the occupation numbers of all fleld modes. 
For the Polyakov loop, only the zero-mode (condensate) 
is considered here. The equation of motion for £ is 
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This equation follows from (|^) when fluctuations of the £ 
fleld are disregarded, as in the present qualitative analy- 
sis. Such an approximation is valid for short times, when 
the energy density is dominated by the zero-mode of the 
Polyakov loop and the contributions from fluctuations of 



£ and of the chiral flelds (the produced mesons) are not 
large. 

The fluctuations of the chiral flelds at any given instant 
are computed from their mode functions, 
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The contribution from the vacuum fluctuations diverges. 
The integral is regularized using the UV cutoff A = 
Nk^k = IO/tt — 1 GeV, which is a sensible value for 
our effective theory. In practice, the dependence on this 
cutoff was found to be weak; if A exceeds all other mass 



scales, the dependence on A is only logarithmic 13 22[|. 

In the above equation, A4 denotes the occupation 
number of a mode with given isospin and with momen- 
tum fc, subtracted for vacuum fluctuations, 



A4 
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and flk = ^JW^+rri^^^^+T^^^. 11^^ is the tadpole self- 
energy of meson field tt^. With the approximations that 
lead to the interaction (0), it is obtained by taking two 
functional derivatives of Cint with respect to tt^, which 
is trivial since >Ci„t is quadratic in the dynamical flelds. 
The expressions for -|- can simply be read off 



from the r.h.s. of eqs. (B|lC|Jll|). 

For the initial condition at t = 0, we consider the case 
that only adiabatic quantum fluctuations of the meson 
flelds exist, but no on-shell particles, i.e. 



(15) 



To calculate the particle production dynamically, we 
make some simplifications. The dynamical picture is that 
of a quench to a temperature T;j, where the condensate 
for £ decays into hadrons. (T/j is not to be confused with 
the decoupling temperature of the chiral field $ from the 
£ field.) The quench approximation seems reasonable as 
the effective potential changes very rapidly around Tc. 
Dynamical simulations of the long wavelength modes on 
a lattice show that the ^-condensate can be quenched to 
perhaps a few percent below Tc [ p^ . Quenching to lower 
temperature is hardly possible, since the £ fleld eventually 
just "rolls down" towards the conflned minimum. This 
is also evident from the rapidly steepening potential in 
Fig. I 

In Fig. ^, we show the result of such a computation, us- 
ing a hadron production temperature of Th — 0.96Tc and 
an initial value for the expectation value of the Polyakov 
loop of £{t = 0) = 0.7. The Afield "rolls down" towards 
the minimum at ^ = 0, and thereby the energy of the de- 
confined state is converted into physical hadronic states. 
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Initially, as few hadrons have been produced, the rela- 
tive number of kaons and pions fluctuates. However, it 
quickly settles to about K/i: sa 0.15. The relative sup- 
pression of kaon production is due to its larger mass. At 
Th = 0.96Tc, the mass for the real part of ^ (defined as the 
curvature of the potential about ^ = 0), « 400 MeV, 
is between the vacuum masses of the pion and the kaon, 
hence the phase space for decay into kaons is smaller than 
for the decay into pions. 

The nearly perfect agreement of our result on the K/tt 
ratio with experimental data from the CERN-SPS |2^] 
and BNL-RHIC accelerators is certainly a numeri- 
cal coincidence. For example, variations of the coupling 
constants and in particular of Th will certainly change 
the result somewhat. Since this is not our major con- 
cern here, we refrain from a more detailed discussion. 
Rather, the main point is that from the behavior of the 
effective Lagrangian for the Polyakov loop depicted in 
Fig. H one is led to expect that the i field decays into 
hadrons at Th very close to T^. At such temperature 
then, < m^{Th) < rriK and so the available phase 
space leads to much fewer kaons than pions. 
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FIG. 3. The K"^ j-K"^ (solid line) and the 77/7r° (dashed 
line) ratios as functions of time. 

To illustrate this fact, we repeat the computation, low- 
ering Th by hand. While this is unrealistic because it 
assumes that the condensate for ^ retains a large expec- 
tation value even below Tc, it shows that as Th is lowered, 
the decay of the condensate for I leaves us with essentially 
equal numbers of kaons and pions. This is due to the in- 
crease of rrif, which makes phase-space restrictions less 
relevant. Fig. ^ confirms this expectation. If the hadrons 
were produced at T = O.STTc, where mi k, 740 MeV is 
above the kaon vacuum mass, we obtain Kjix k, 0.8. 
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FIG. 4. The K"^ I-k"^ (solid line) and the ry/Tr" (dashed 
line) ratios as functions of time when T/Tc = 0.87. 

The theoretical expectation that the decay of the con- 
densate for i into hadrons occurs close to Tci which was 
based on the rapid variation of the effective potential for 
I from Fig. I p^Jl^ , is thus perfectly consistent with the 
constraints from experimental data on the Kjix ratio in 
high-energy heavy-ion scattering at the CERN-SPS |2^] 
and BNL-RHIC . What would perhaps be even more 
interesting from the point of view of learning about the 
presence of fundamental QCD scales in the hadroniza- 
tion process is to look for situations where our calcula- 
tion breaks down. Fortunately, there indeed seems to be 
an experimentally accessible regime where this can be 
tested, namely high-energy, high- multiplicity pp inelastic 
reactions. 

At given energy y^/A, and at central rapidity, the K/t: 
ratio in heavy-ion collisions approaches that from pp (or 
pp) reactions at the same ^/s as the energy increases . 
A natural question to ask is whether indeed they asymp- 
totically approach the same limit. In fact, one may even 
compare central heavy- ion collisions at, say, BNL-RHIC 
energy to high-multiplicity pp (or pp) at much higher en- 
ergies, where the charged multiplicity per participant and 
the energy density at central rapidity are comparable to 
those from heavy-ion collisions, and where one may very 
well produce a transient QGP state. In the absence of a 
length scale l/mg, one would naively expect that eventu- 
ally, at sufficiently high energy, the K/t: ratio from high- 
multiplicity pp approaches that from central heavy-ion 
collisions at BNL-RHIC energy (also if high-multiplicity 
fluctuations in pp are due to jets). 

On the other hand, if the free energy of QCD around 
Tc is dominated by a condensate for the Polyakov loop, 
hadron production is governed by a scale mi as explained 
above. Thus, hadronization from the decay of the con- 
densate is associated with a length scale 1/mi. It ap- 
pears quite reasonable to assume that it can only dom- 
inate hadron production in systems of size greater than 
1/mi. (In small systems, on the other hand, fluctua- 
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tions dominate.) With rrii somewhere between and 
rriK, as given by the potential for £ from Fig. ||, this con- 
dition may not hold for pp reactions, regardless of how 
high an energy density is achieved. Note that we are not 
referring here to thermodynamical effects like for exam- 
ple canonical suppression due to exact conservation of 
strangeness (which becomes less important for very high 
energy and high-multiplicity pp, when there is more than 
one charged kaon pair per unit of rapidity). We are re- 
ferring to whether or not one can probe, in principle, the 
deconfined phase of the non-abelian SU{3) gauge the- 
ory in a system of linear dimension l/nii (or less), at a 
temperature T T^- 

Consider pure gauge theory with three colors, for which 
good lattice data exists. In d Euclidean dimensions, the 
two-point function for the real part of the Polyakov loop 
behaves as 



{lr[r%{0)) - {irf / d'^k 



/TO,XW2-1) 

~ j ■ (16) 

The correlation function for the imaginary part of £ 
falls off more rapidly, corresponding to a larger mass for 
li For d 3, Ki/2{mir) ^ I / y/mPr exp{mi r) , 

hence 



{£Ar)£rm - {£r 



(17) 



When the size of the system r < ^/rni, the correlation 
function ( |l7|) exhibits a power-law behavior 



{£r{r)£rm ~ {£r) 



(18) 



for all "accessible" distances. The r— dependence of the 
subtracted correlation function ( p^ ) is in fact that for 
an abelian pure gauge theory (photons), namely just the 
Coulomb potential. The fact that the Polyakov loop 
aquires a mass m^, which is the feature of the non- 
abelian gauge theory, can only be probed on distance 
scales r ^ l/mg. Our computation of hadron production 
applies to systems larger than I /mi, when the line-line 
correlation function falls off exponentially and when the 
free energy is dominated by the condensate for £. Qual- 
itatively, hadron production about Tc in systems of size 
1/me might differ from our model. 

Data on the K/n ratio at central rapidity from pp 
reactions at ^/s — 1800 GeV at the Fermilab-Tcvatron 
have been analysed as a function of charged multiplic- 
ity 1 25 2^. For the highest-multiplicity events, corre- 
sponding to dNch/drj w 25 the initial energy density at 
proper time tq = 1 fm/c was estimated to be 6 GeV/fm'^, 
comparable to the average value from central Au-|-Au col- 
lisions at BNL-RHIC energy y/^/A ^ 100-200 GeV 
A first analysis p^] seemed to indicate that K/ n increases 



as a function of event multiplicity, from about 10% to 
about 14%, nearly the value from Au-|-Au at RHIC. How- 
ever, error bars for the most inelastic events were large, 
and a second analysis [ p6[ did not confirm the rising K/t: 
ratio but showed it to be completely flat as a function of 
dNch/drj. A new analysis of this and other observables 
with the colliding proton beams at RHIC would be most 
interesting. 

In summary, we have discussed dynamical hadron pro- 
duction at the confinement transition, which is realized 
here as the restoration of the global Z(3) symmetry for 
the Polyakov loop £. Specifically, the proposal to under- 
stand QCD near in terms of a condensate for £ 13 
predicts that hadron production occurs very near T^. 
This follows from the rapid change of the effective po- 
tential for the Polyakov loop about which we deduce 
from lattice data on the finite-temperature string ten- 
sion, pressure, and energy density. From our effective 
Lagrangian, the mass of £ just below Tc is greater than 
the vacuum mass of the pion but less than that of the 
kaon. Therefore, qualitatively one expects about an or- 
der of magnitude more pions than kaons to be produced, 
as seen in experimental data from CERN-SPS and BNL- 
RHIC. 

The hadronization process is characterized by a length 
scale 1/m^, which is on the order of 1/2 — 1 fm, just 
below Tc- Thus, in small systems hadronization may 
not occur from the decay of a condensate for £, regard- 
less of the initial energy density. Differences in the K/tt 
ratio (and in other relative hadron multiplicities) from 
high-multiplicity pp to heavy-ion collisions at the SPS 
or RHIC [||,^ could perhaps be due to a condensate 
scale 1/mi. If such a point of view can be established 
experimentally, it could provide further support for the 
formation of a deconfined phase in collisions of heavy ions 
at high energies. 
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